Damage threshold of LiF crystal irradiated by femtosecond hard XFEL pulse sequence.

Here we demonstrate the results of investigating the damage threshold of a LiF crystal after irradiating it with a sequence of coherent femtosecond pulses using the European X-ray Free Electron Laser (EuXFEL). The laser fluxes on the crystal surface varied in the range ∼ 0.015-13 kJ/cm2 per pulse when irradiated with a sequence of 1-100 pulses (tpulse ∼ 20 fs, Eph = 9 keV). Analysis of the surface of the irradiated crystal using different reading systems allowed the damage areas and the topology of the craters formed to be accurately determined. It was found that the ablation threshold decreases with increasing number of X-ray pulses, while the depth of the formed craters increases non-linearly and reaches several hundred nanometers. The obtained results have been compared with data already available in the literature for nano- and picosecond pulses from lasers in the soft X-ray/VUV and optical ranges. A failure model of lithium fluoride is developed and verified with simulation of material damage under single-pulse irradiation. The obtained damage threshold is in reasonably good agreement with the experimentally measured one.

Mass and heat transfer between evaporation and condensation surfaces: Atomistic simulation and solution of Boltzmann kinetic equation.

Boundary conditions required for numerical solution of the Boltzmann kinetic equation (BKE) for mass/heat transfer between evaporation and condensation surfaces are analyzed by comparison of BKE results with molecular dynamics (MD) simulations. Lennard-Jones potential with parameters corresponding to solid argon is used to simulate evaporation from the hot side, nonequilibrium vapor flow with a Knudsen number of about 0.02, and condensation on the cold side of the condensed phase. The equilibrium density of vapor obtained in MD simulation of phase coexistence is used in BKE calculations for consistency of BKE results with MD data. The collision cross-section is also adjusted to provide a thermal flux in vapor identical to that in MD. Our MD simulations of evaporation toward a nonreflective absorbing boundary show that the velocity distribution function (VDF) of evaporated atoms has the nearly semi-Maxwellian shape because the binding energy of atoms evaporated from the interphase layer between bulk phase and vapor is much smaller than the cohesive energy in the condensed phase. Indeed, the calculated temperature and density profiles within the interphase layer indicate that the averaged kinetic energy of atoms remains near-constant with decreasing density almost until the interphase edge. Using consistent BKE and MD methods, the profiles of gas density, mass velocity, and temperatures together with VDFs in a gap of many mean free paths between the evaporation and condensation surfaces are obtained and compared. We demonstrate that the best fit of BKE results with MD simulations can be achieved with the evaporation and condensation coefficients both close to unity.

Dynamic fracture of tantalum under extreme tensile stress.

The understanding of fracture phenomena of a material at extremely high strain rates is a key issue for a wide variety of scientific research ranging from applied science and technological developments to fundamental science such as laser-matter interaction and geology. Despite its interest, its study relies on a fine multiscale description, in between the atomic scale and macroscopic processes, so far only achievable by large-scale atomic simulations. Direct ultrafast real-time monitoring of dynamic fracture (spallation) at the atomic lattice scale with picosecond time resolution was beyond the reach of experimental techniques. We show that the coupling between a high-power optical laser pump pulse and a femtosecond x-ray probe pulse generated by an x-ray free electron laser allows detection of the lattice dynamics in a tantalum foil at an ultrahigh strain rate of [Formula: see text] ~2 × 108 to 3.5 × 108 s-1. A maximal density drop of 8 to 10%, associated with the onset of spallation at a spall strength of ~17 GPa, was directly measured using x-ray diffraction. The experimental results of density evolution agree well with large-scale atomistic simulations of shock wave propagation and fracture of the sample. Our experimental technique opens a new pathway to the investigation of ultrahigh strain-rate phenomena in materials at the atomic scale, including high-speed crack dynamics and stress-induced solid-solid phase transitions.

Solitary Nanostructures Produced by Ultrashort Laser Pulse.

Laser-produced surface nanostructures show considerable promise for many applications while fundamental questions concerning the corresponding mechanisms of structuring are still debated. Here, we present a simple physical model describing those mechanisms happened in a thin metal film on dielectric substrate irradiated by a tightly focused ultrashort laser pulse. The main ingredients included into the model are (i) the film-substrate hydrodynamic interaction, melting and separation of the film from substrate with velocity increasing with increase of absorbed fluence; (ii) the capillary forces decelerating expansion of the expanding flying film; and (iii) rapid freezing into a solid state if the rate of solidification is comparable or larger than hydrodynamic velocities. The developed model and performed simulations explain appearance of microbump inside the focal spot on the film surface. The model follows experimental findings about gradual transformation of the bump from small parabolic to a conical shape and to the bump with a jet on its tip with increasing fluence. Disruption of the bump as a result of thinning down the liquid film to a few interatomic distances or due to mechanical break-off of solid film is described together with the jetting and formation of one or many droplets. Developed theory opens door for optimizing laser parameters for intended nanostructuring in applications.

Laminar, cellular, transverse, and multiheaded pulsating detonations in condensed phase energetic materials from molecular dynamics simulations.

The development of condensed-phase detonation instabilities is simulated using moving window molecular dynamics and a generic AB model of a high explosive. It is found that an initially planar detonation front with one-dimensional flow can become unstable through development of transverse perturbations resulting in highly inhomogeneous and complex two- and three-dimensional distributions of pressure and other variables within the detonation front. Chemical reactions are initiated in localized transverse shock fronts and Mach stems with a pressure and temperature higher than those predicted by classic Zel'dovich, von Neumann, and Doering detonation theory. The two-dimensional cellular and transverse and three-dimensional pulsating detonation structures are found by varying the physico-chemical properties of AB energetic material, sample geometry, and boundary conditions. The different regimes of condensed-phase detonation that can develop from instabilities within a planar detonation front exhibit structures, although at a much smaller scale, that are similar to those observed in gases and diluted liquids.

Evolution of shock-induced orientation-dependent metastable states in crystalline aluminum.

The evolution of orientation-dependent metastable states during shock-induced solid-liquid phase transitions in crystalline Al is followed using moving window molecular dynamics simulations. The orientation-dependent transition pathways towards an orientation-independent final state Hugoniot include both "cold melting" followed by recrystallization in [110]- and [111]-oriented shock waves and crystal overheating followed by melting in [100] shock waves. The orientation-dependent dynamics take place within a zone that can extend up to hundreds of nanometers behind the shock front.

Two-zone elastic-plastic single shock waves in solids.

By decoupling time and length scales in moving window molecular dynamics shock-wave simulations, a new regime of shock-wave propagation is uncovered characterized by a two-zone elastic-plastic shock-wave structure consisting of a leading elastic front followed by a plastic front, both moving with the same average speed and having a fixed net thickness that can extend to microns. The material in the elastic zone is in a metastable state that supports a pressure that can substantially exceed the critical pressure characteristic of the onset of the well-known split-elastic-plastic, two-wave propagation. The two-zone elastic-plastic wave is a general phenomenon observed in simulations of a broad class of crystalline materials and is within the reach of current experimental techniques.